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^_^ ' Abstract 

i-^ ' We study the prospects of measuring the CKM matrix element |Vfs| at the LHC with the top quarks 

Ph' _ _ __ 

' I produced in the processes pp — > ttX and pp — > t/tX, and the subsequent decays t -^ W^s and t ^ W s. 

D I To reduce the jet activity in top quark decays, we insist on tagging the W^ leptonicahy, W^ — )■ i^vg 

1—^, {£ — e,^,,T), and analyse the anticipated jet prohles in the signal process t — > Ws and the dominant 

r-fs , background from the decay t — ?> Wh. To that end, we analyse the Vi) {K'~' and A) distributions in the s- and 

jj^ . 6-quark jets concentrating on the energy and transverse momentum distributions of these particles. The 

^^ I l^Os emanating from the t — >■ Wb branch have displaced decay vertexes from the interaction point due to the 

■^^ ■ weak decays b ^)- c ^)- s and the 6-quark jets are rich in charged leptons. Hence, the absence of secondary 

l/^ ' vertexes and of the energetic charged leptons in the jet provide additional (6-jet vs. s-jet) discrimination in 

^— s ' top quark decays. These distributions are used to train a boosted decision tree (BDT), a technique used 

^~~^ ' successfully in measuring the CKM matrix element |Vtf,| in single top production at the Tevatron. Using 

^^ ' the BDT classifier, and a variant of it called BDTD, which makes use of decorrelated variables, we calculate 

KS ' the BDT(D)-response functions corresponding to the signal (t -^ Ws) and background (t -^ Wb). Detailed 

Cu ^ simulations undertaken by us with the Monte Carlo generator PYTHIA are used to estimate the background 

rejection versus signal efficiency for three representative LHC energies yi = 7 TeV, 10 TeV and 14 TeV, 

of which only the analysis for the y/s — 14 TeV case is shown in detail. We argue that a benchmark with 

10% accuracy for the signal {t —>■ Ws) at a background (i — > Wb ) rejection by a factor lO'^ (required due 

to the anticipated value of the ratio l^tsp/lVibp ~ 1.6 x 10^^) can be achieved at the LHC@14 TeV with 

an integrated luminosity of 10 fb^^. 
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I. INTRODUCTION 

It is now fifteen years that the top quark was discovered in proton-antiproton colhsions at 
the Tevatron [l|, \^. Since then, a lot of precise measurements have been undertaken at the two 
Fermilab experiments, CDF and DO. Among the highhghts are the measurements of the top quark 
mass, currently having an accuracy of about 0.75%, the ti production cross section with about 9% 
accuracy [3], and the observation of the electroweak single top production \4,\^- Of these, the single 
top (or anti-top) production cross section depends on the charged current couplings tqW, where 
q = d,s,b, which in the standard model (SM) are governed by the Cabibbo-Kobayashi-Maskawa 
(CKM) quark mixing matrix Vckm [a, l3]: 

Vud ^us Kb' 

Vckm = \ Vcd Vcs Vcb 
Vtd Vts VtJ 

In the Wolfenstein Parametrisation [8[, this matrix is expressed as 

/ l-iA^ A A\^p-iri)\ 



Vckm 



-X{l + iA^X'^r]) 1-^A2 AA2 

\A\^ (1 - P - ii]) -AX^ (l + iX^r]) 1 / 



where A, A, p and t] are the Wolfenstein parameters. 

The cross section a{pp — )• t/iX) has provided the first direct measurement of the dominant 
CKM-matrix element \Vtb\- In this analysis, it is assumed that the CKM matrix elements \Vtd\ and 
I Vts I are much smaller than | Vtb \ , but no assumption is made about the unitarity of the 3x3 CKM 
matrix. To obtain |Vtfcp, the measured cross section for an suumed top quark mass is divided by 
the theoretical cross section for \Vth\ = 1- Following this procedure, the CDF measurements yield 
\Vtb\ = 0.91±0.11(stat +syst)±0.07(theory), which in turn gives \Vtb\ > 0.71 at 95% C.L. p. The 
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corresponding limit from DO is \Vtb\ > 0.78 at 95% C.L. J4|. The combined CDF and DO analysis 
assumes 170 GeV as the top quark mass and yields \Vtb\ = 0.91 ± 0.08 with \Vtb\ > 0.79 at 95% 
C.L. using a{pp -^ t/i+X) = 3.14 pb [lot. A theoretical value a{pp -^ t/t + X) = 3.46 pb [g| as 
input yields \Vtb\ = 0.88 ± 0.07 with \Vtb\ > 0.77 at 95% C.L. [11]. There also exist limits on this 
matrix element obtained from the decays of the top quarks by tagging the b-quark jet in the final 
state. Defining the ratio Rt^quiv j^,^_^^yy^ ^\'^^^l ^^^^^f^y^^K = iVjfep, where use has been made of 
the CKM unitarity in the second equality, CDF and DO measurements yield \Vtb\ > 0.78 [12] and 
iVtfel > 0.89 [131, respectively. 

The above determination of the matrix element \Vtb\, obtained from the direct single top produc- 
tion and the 6-tagged decays of the top quark, can be compared with the indirect determination of 
the same based on a number of loop-induced processes in which top quark participates as a virtual 



state, such as the B^- B^ and B'^- B^ mixings, the radiative decay B — t- Xg^ and the CP-violation 
parameter ex in the Kaon sector. Overall fits of the CKM unitarity yield, comparatively speaking, 
an infinitely more accurate value \Vtb\ = 0.999133(44) ij]. This precision, in all likelihood, will 
not be matched by the direct determination of \Vth\, as experiments at the LHC are expected to 
reach an accuracy of a few per cent on this quantity. Neverthelss, a determination of \Vtb\ with 
such an accuracy will be very valuable to constrain beyond-the-SM physics models. A good case 
in point is a model with four generations, in which \Vtb\ can be as low as 0.93 |15l ]. 

We go a step further and explore in this paper the prospects of measuring the matrix element 
\Vts\ at the LHC. In the Wolfenstein parametrisation [S], this matrix element is given as \Vts\ = 
AX^ + 0{X'^). The best-fit values from the unitarity fits are: A = 0.814, A = 0.2257, yielding \Vts\ = 
0.0407 lb 0.001 iJ]. The smallest matrix element in the third row of the CKM matrix is Vtd, and 
its value from the CKM unitarity fit is posted as 114^1 = AX^y^{l - pY + rf = (8.74J:o:37) x lO'^. 
Direct determination of these matrix elements will require a good tagging of the i — )■ s transition 
(for \Vts\) and t ^ d transition (for \Vtd\) in the top quark decays, and a very large top quark 
statistics, which will be available only at the LHC in the foreseeable future from the processes 
pp ^ tt + X and pp — )• t/t + X. Just as for the direct measurement of \Vtb\, there is also a lot of 
interest in the direct measurements of Vts and Vtd, as the absolute values of these CKM matrix 
elements can be modified by approximately a factor 2 from their SM values quoted above, taking 
the example of a four-generation extension of the SM [15]. Lacking a good tagging for the t — >• d 
transition, and also because of the small size of the CKM-matrix element, \Vtd\ = 0(10^^), we 
concentrate here on the direct measurements of \Vts\ at the LHC. 

In order to be able to measure \Vts\ directly, one has to develop efficient discriminants to 
suppress the dominant decay t ^ W b. As the first step, we propose to tag only those events in 
which the W^ decay leptonically to reduce the jet activity in top quark decays. The emerging 
s-quark from the top quark decay t ^ W s, and the collinear gluons which are present in the 
fragmentation process anyway, will form a hadron jet. We suggest tagging on the yOs {K^s and 
As) in this jet, and measure their energy and transverse momentum distributions. Energetic yOs 
are also present in the 6-quark jets initiated by the decay t ^ W b and the subsequent weak decays 
6 — 7- c — >• s. However, in this case, the T^Os will be softer, will have displaced vertexes (from the 
interaction point) and they will be often accompanied with energetic charged leptons due to the 
decays b — )■ £^X. Absence of a secondary vertex and paucity of the energetic charged leptons in 
the jet provide a strong discrimination on the decays t — )■ Wb without essentially compromising 
the decays t ^ W s. Thus, the scaled energy and transverse momentum distributions of the K^s, 
As and is, and the secondary decay vertex distributions (dN/dr) are the quantities of principal 
interest. Here r is the distance traversed in the transverse plane, i.e. the plane perpendicular to 



the beam axis or r — cf) plane, by the 6-quark before decaying, smeared with a Gaussian resolution 
to take into account realistic experimental conditions. 

We have assumed two representative r.m.s. values (cr(vertex) =1 mm and 2 mm) for the Gaus- 
sian, where 2 mm is more realistic. Experimentally, b-tagging algorithms are based on measure- 
ments of the impact parameter from the B meson charged tracks. The power separation between 
b- and light-jets using this so called 2D-method is similar to the one we used in our analysis with 2 
mm resolution. We also show results with 1 mm resolution in order to illustrate how important it 
will be for the Super LHC (SLHC) to improve on the b-tagging efficiency. These distributions are 
calculated for the processes pp — )• tiX and pp — )• t/iX, for the signal (t — )• VF s) and background 
{t^W b). 

Having generated these distributions, characterising the signal t ^ W s and the background 
t ^ W b events, we use a technique called the Boosted Decision Tree (BDT) - a classification 

n 

model used widely in data mining [16f] ~ to develop an identifier optimised for the t ^- W s 
decays. In our calculation, we use both BDT and a variant of it called BDTD (here D stands for 
decorrelated) , where possible correlations in the input variables are removed by a proper rotation 
obtained from the decomposition of the square root of the covariance matrix, to discriminate the 
signal events from the large backgrounds. We recall that this technique has been successfully used 
to establish the single top quark production in pp collisions at the Tevatron [J, |5( (see [l7| for 
details). Briefiy, the generated input is used for the purpose of training and testing the samples. 
We provide the input in terms of the variables discussed earlier for the signal (t — )• Ws) and the 
background (t — t- Wb), obtained with the help of a Monte Carlo generator. This information is 
used to develop the splitting criteria to determine the best partitions of the data into signal and 
background to build up a decision tree (DT). The separation algorithm used in splitting the group 
of events in building up DT plays an important role in the performance. The software called the 



Toolkit for Multivariate Data Analysis in ROOT (TMVA) [l8|] is used for the BDT(D) responses 
in our analysis. Detailed simulations presented here are done using PYTHIA [19(] to model the 
production processes, gluon radiation, fragmentation and decay chains, and the underlying events. 
We calculate the signal {t ^ W s) efficiencies for two cases called bb/bs and bs/ss (defined in 
section II) for an assumed (Gausian) vertex smearing with an r.m.s. value of 2mm and 1mm. 
Concentrating on the bb/bs case, when only one of the top (or antitop) quark decays via t — )• sVF"*", 
compared to the case when both the top and antitop decay via the dominant transition t — )• bW~^, 
these efficiencies lie typically between 5% (for the 2mm smearing) and 20% (for the 1mm case) for 
a background (t ^- W b) rejection by a factor 10'^ (see Table I). 

Note that this level of background rejection is necessary due to the anticipated value of the 
ratio iVfsp/lVihp ~ 1.6 x 10^^. The required integrated LHC luminosity to determine \Vts\ directly 



is estimated as 10 fb~^ at 14 TeV. Numerical analysis reported here is carried out for three repre- 
sentative LHC energies: ^/s = 7 TeV, 10 TeV and 14 TeV, but we present the detailed results only 
for ^/s = 14 TeV as the distributions for y/s = 7 TeV and 10 TeV are similar to the 14 TeV case. 

In section 2, we study the process pp — )■ ttX, reviewing first the production cross sections 
at the LHC energies. The energy-momentum profiles of the signal {t ^ W s) and background 
{t — ;■ Wh) events produced in the tt pair production process and the subsequent decays pp — )■ 
t{-^ W+b) i{-^ W-b)X, pp -^ t{-^ W+b) t{^ W-s)X, pp -^ t{^ W+s) t{^ W'b)X and 
pp — )• t(— )• VF~'~s) t(— )• W~s)X are worked out. The last of these has a very small branching ratio 
and its measurement would require a huge LHC luminosity (we included this case for the sake of 
completeness). Tagging efficiencies for pp — )• tiX calculated with the BDT(D) classifier are shown 
in Table 1. Numerical results in this table are presented as (bb/bs) and (bs/ss), corresponding to 
the cases when only one of the t (or t) decays via t — )• W'^ s (or t — )• W~ s) and when both the t 
and t decay via t — )• W~^ s and i — )• W~ s, respectively. 

Section 3 is a repeat of the above analysis for the single top production process pp — )• t/tX at 
the LHC. The end-product of this analysis chain is again the background rejection vs. the signal 
efficiency based on the BDT(D) response functions. The tagging efficiencies, calculated for ^/s = 7, 

10, and 14TeV with the BDTD classifier are presented in Table 2. Section 4 briefly summarises 
our results and outlook. 

11. ANALYSIS OF THE PROCESS pp -^ ttX AND THE SUBSEQUENT DECAYS t -^ 

Wh,Ws 

Theoretical predictions of the top quark production at the LHC have been obtained by including 
up to the next-to-next-to- leading order (NNLO) corrections in the strong coupling constant 22I-25I. 
They have been updated taking into account modern parton distribution functions (PDFs) 



A typical estimate is: cr(pp -^ ttX) = 874+33pb for nit = 173 GeV and ^/s = lA TeV [28(], where 
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the errors reflect the combined uncertainties in the factorisation and normalisation scales and in 



the parton distribution functions (PDF). Other independent NNLO calculations yie! 



d similar cross 



section, though the error budgeting is somewhat different. Kidonakis and Vogt 2M P^^ the cross 
section a{pp — )• ttX) = 894 it 4(kinematics)t44(scale)l3;^(PDF) for the same values of mt and ^/s, 
using the CTEQ6.6M PDFs [271, and (T{pp -^ tiX) = 943 ± 4(kinematics)+^^(scale) ± 12(PDF), 
using the MRST 2006 PDFs [26J]. Compared to the tt production cross section at the Tevatron 
{^/s = 1.96 TeV), (j{pp — )■ tiX) = 7.34lQ'3g pb [28], one expects a rise in the tt cross section by more 
than two orders of magnitude between the Tevatron and the LHC@14 TeV. The cross sections at 
the lower LHC energies, 7 and 10 TeV, have also been calculated 
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28(], with a{pp -^ ttX) ~ 400 



pb at 10 TeV and about half that number at 7 TeV. Thus, for the top quark physics, the dividends 
in going from 7 to 14 TeV are higher by a good factor 4. 

For the numerical results shown here we have used the PYTHIA Monte Carlo [l9[ to generate 10^ 
events for the process pp — t- ttX, followed by the decay chains t — t- W~^h, W'^s and i — )• W~b, W~s. 
As stated in the introduction, the W^ are forced to decay only leptonically W^ — )• l^fi {£ = e, //, r) 
to reduce the jet activity from the non-leptonic decays of the W . This corresponds to an integrated 
luminosity of 10 fb^^ at 14 TeV. For an estimated efficiency of 5% at a 10~^ background rejection, 
and |14sP ^ 1-7 x 10^^, as in the SM, this means that we expect 0.05 x 2 x 1.7 x 10"^ x 10*^ = 170 
signal events with a background of 10^^ x 10^ = 10^ events, giving a significance of 170/\/l000 i.e. 
more than 5a. 

We then concentrate on the VO production, which for the experimental conditions at the two 



main detectors ATLAS [20'1 and the CMS [2l|] imphes VO = K'^ or VO = A, as the long-lived Kl 



will decay mostly out of the detectors. However, both K^ and A can be detected by ATLAS and 
CMS and their energy and momentum measured with reasonably good precision. In the present 



analysis, we reconstruct FOs and soft leptons in the rapidity range |r/| < 2.5 20]. In addition, 
we require the yO's decay radius to lie in the range 20 to 600 mm. These acceptance cuts are 
acceptable for both multipurpose detectors mentioned above, and they will be used in the analysis 
described in this and the next section. 

We will show the distributions for ^/s = 14 TeV, the designed LHC center-of-mass energy. The 
X'^-energy distribution is shown in the left-hand frame in Fig [T] plotted as a function of the scaled 
energy Xk = Ex/Ejct- For this study, the jet energy is set equal to the quark energy produced 
in the decay t — )• Wb, Ws. In a realistic simulation of the experimental measurements, one would 
require a functional definition of the jet, for example using an angular cone, which will then define 
the jet energy -Ejet, and hence xk- The transverse momentum of the K^s, pt{K^) (in GeV), is 
shown in the right hand frame in Fig[TJ In both of these frames, the solid histograms correspond 
to the decay t ^ W s and the dashed ones to the decay t ^ W b. As expected, the decay chain 
t — )• W s{-^ K^) has a much stiffer distribution both in Xk and pt{K^), as the K^^s descending 
from the decay chain t ^ W 6(— )• c — )• s) are rapidly degraded in these variables due to the 
subsequent weak decays. The corresponding distributions for the As are shown in the lower two 
frames in Fig. [TJ They are qualitatively very similar to those of the K^s. 

We now show the distributions in the charged lepton energy from the decays t — )• 6 — )• ^^A" 
and t — )• s — )• i^X in Fig. [2l showing the scaled lepton energy in the variable X^ = Ei/Ej^t 
(upper frame) and in p^, the transverse momentum of the charged leptons (middle frame). This 
distribution quantifies the richness of the 5-jets in charged leptons and the stiff character of the 
energy /transverse momentum distributions due to the weak decays, as compared to the leptons 
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FIG. 1: pp — > ttX at ^/s ~ 14 TeV. Upper left frame: scaled-i^"-energy distributions dN/dxK from 
t ^f W s{-^ K^X) (solid histogram) and t —>■ W 6(— > K'^X) (dashed histogram). Upper right frame: 
Transverse momentum distributions of the K'^s measured w.r.t. beam axis dN/dpT^ (in GeV) in the same 
production and decay processes as in the left frame. Lower frames show the distributions dN/dx\ and 
dN/dpT^ (in GeV) for t ^ W s{-^ AX) (solid histogram) and t^W b{-^ AX) (dashed histogram). 



from s — )• £ X, which are all soft and coming from the leptonic decays of the various resonances 
produced in the fragmentation of the s-quark. Absence of energetic charged leptons in the s-quark 
jet in the decay t — )• VF s is a powerful tool in reducing the background from the otherwise much 
more prolific process t ^ W b. The final set of distributions from our Monte Carlo simulation is 
the secondary decay vertex distribution (lower frame) , smeared with a Gaussian distribution with a 
r.m.s. of 2 millimetres, shown in terms of a variable called r (measured in millimetres). The decay 
length for the t ^ W b case is calculated as 'Jcti,, where 7 is the Lorentz factor, and ctj, = 0.45 
mm, corresponding to an average 6-quark lifetime taken as tj, = 1.5 ps from the PDG fl4|. This 
distribution, which reflects the long lifetime of the 6-quark (respectively of the B and A^ hadrons), 
as opposed to the lack of a secondary vertex from the s-quark fragmentation process, is also a very 
powerful discriminant of t ^ W b vs. t ^ W s decays. 

Having generated these distributions, characterising the signal t ^ W s and background t — )• 
W b events in the process pp — )■ tiX at the LHC, we use the BDT and BDTD classifiers, discussed 
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FIG. 2: pp -^ ttX at ^/s = 14 TeV. Upper frame: scaled-i?^-energy distributions, dN/dxi, from t ^- W s{-^ 
i^X) (solid histogram) and t ^ W b{-^ i^X) (dashed histogram). Middle frame: Transverse momentum 
distributions of the ^ s measured w.r.t. beam axis, dN/dpT^ (in GeV), in the same production and decay 
processes as in the upper frame. Lower frame: Secondary decay vertex distributions in the variable r (in 
millimetres) for the two decay chains t -^ W s (solid histogram) and t ^ W b (dashed histogram), obtained 
by smearing the decay length with a Gaussian having an r.m.s. value of 2 mm. 



in the introduction. In Fig. [3] (left frame), we show the BDTD response functions, showing that a 
clear separation between the signal {t ^ W s) and background {t ^ W b) events has been achieved. 
The background rejection vs. signal efficiency for the pp — )• ti events is shown in Fig. [3] (right frame) 
for both the BDT and BDTD classifiers, which give very similar results. The evaluation results 
ranked by the best signal efficiency and purity are shown numerically in Table 1. The entries in 
this table show that a background rejection of 10'^ can be achieved at a signal efficiency of about 
5% to reach the SM-sensitivity of the CKM matrix element \Vts\- The statistical uncertainty in 
this efficiency is 0.026%. It goes up to 0.067% for efficiency values at the upper end. 
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FIG. 3: pp -^ ttX at y/s — 14 TeV. Left frame: The normalised BDTD response, calculated by using the 
TMVA (see text). The signal (dark shaded) from the decay t ^ W s and the background (light shaded 
with dotted lines) from the decay t ^ W b are clearly separated in this variable. Right frame: Background 
rejection vs. signal efficiency calculated from the BDT(D) response. The result using the BDT classifier is 
very similar and hardly distinguishable from the one obtained with the BDTD response. 



The distributions at ^/s = 7 and 10 TeV are very similar to the corresponding ones shown 
in Fig. [1] for 14 TeV. Hence, the characteristic differences that we have shown at ^/s = 14 TeV 
emanating from the top quark decays t — )• Wb and t — )• Ws in the VO and charged lepton energy- 
and transverse momentum spectra are also present at the lower energies. 

Based on the above analysis we have calculated the tagging efficiencies for the decay t ^- W s 
(signal) for an acceptance of 0.1% for the decay t — )• Wb (background). The acceptance level is 
motivated by the anticipated value of the ratio of the t — t- Ws and t — )■ Wb decay rates, which in 
the SM is O(10~^). The tagging efficiencies for the three centre-of-mass energies at the LHC (7, 
10 and 14 TeV) are given in Table 1 for two different vertex smearing (1 mm and 2 mm), assuming 
a Gaussian distribution. The entries shown as bb/bs correspond to the comparison for top pair 
production process pp — )• ttX with both the t and t decaying via the dominant process t — t- W~^b 
and i — )• W~b, respectively, and in which only one of the t or i quarks decays via t — )■ W~^s or 



TABLE I: Tagging efficiencies (in %) for the process pp —>■ ttX, foUowed by the decay t —> Ws (signal) and 
t — >■ Wb{ background), calculated for an acceptance of 0.f% for the background at three LHC centre-of- 
niass energies. Two Gaussian vertex smearing (having an r.m.s. values of 2 mm and 1 mm) are assumed for 
calculating the displaced vertex distributions dN/dr. The cases bb/bs and bs/ss are explained in the text. 
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t —7- W~s, and the other decays via t — )• W^b or i — t- W~b (signal events). The entries marked 
as bs/ss correspond to the cases where either the t or t quarks decays via t — )• W~^s or t — )• W~s 
and both t and t quarks decay via t — )• W~^s and t — t- W~s. The branching ratio for the case 
{t — )• W~^s)){i ^ W~s) is exceedingly small, O(10~^). The entries in Table 1 for this case (bs/ss) 
show that at the considerable price of the reduced sample, one can get much better efficiencies. 



III. ANALYSIS OF THE PROCESS pp -^ t/tX AND THE SUBSEQUENT DECAYS t 

Wb,Ws 



The single top (or anti-top) cross sections in hadron hadron collisions have been calculated 



in the NLO approximation 



Id, 



3Cll432l|. Recalling that there are three basic processes at the 
leading order which contribute to CT{pp — )• t/tX), namely the t-channel: qb — )• q't, the s-channel: 
qq' — )• bt; and the associated tW production bg — )• tW~ , the cross sections estimated at the Tevatron 
are [J: a^-^^^^^^^ (pp -^ tX) = a{pp -^ tX) = 1.14 ± 0.06 pb, cr^^'=^^°'^'^'(pp -^ tX) = a{pp -^ 
iX) = 0.53±0.02 pb, and cr^w -c\v!.nne\ ^^^ _^ ^x) = a{pp -^ iX) = 0.14±0.03 pb, putting the single 
top (or antitop) cross section o"(pp — t- tX) = a{pp — t- tX) ~ 1.8 pb at the Tevatron. These cross 
sections have to be compared with the theoretically projected cross sections at the LHC@14 TeV: 
^t-channei^^p ^ ^X) = 149 ± 6 pb, a^-^^^'^^''\pp -^ tX) = 91 ±4 pb, a"-^^'^™^'(pp -^ tX) = 7.71^;^ 
pb, CT^-'=^^°'^'='(pp -^ iX) = 4.3 ± 0.2 pb, and ^jiMz-channei^^p ^ ^x) = a{pp ^ iX) = A3 ± 5 pb. 
Thus, one expects a{pp — t- tX) ~ 200 pb and about half this number for a{pp — t- tX), yielding 
the summed single top and antitop cross sections at about 300 pb at the LHC@14 TeV, also 
approximately two orders of magnitude larger than those at the Tevatron. With a luminosity of 10 
fb^^, one anticipates 0(3 x 10^) single top (or anti-top) events, i.e. O(IO^) events in the leptonic 
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FIG. 4: pp — )■ t/iX at ^/s — 14 TeV. Upper frames: Scaled energy distributions dN/dxK and the 
transverse momentuni distribution dN/dpTK from the decays t — ^ VFs(— > K'^X) (soHd histograms) and 
t — !• Wh{~^ K^X) (dashed histograms). Lower frames: Scaled energy distributions dN/dxA and the 
transverse momentum distribution dN/dpTf^ from the decays t — >■ Ws{-^ AX) (solid histograms) and 
t -> Wb{~^ AX) (dashed histograms). 

channel. Thus, the rise in the cross sections for a single top (or antitop) production between the 
Tevatron and the LHC@14 TeV is also very marked. 

In the analysis shown here, we have again resorted to the Monte Carlo generator PYTHIA, 
which models so far only the s-channel single top production process pp — t- VF — )■ tb. As discussed 
above, this has the smallest (of the three channel) cross section. Moreover, the tW channel process 
bg —7- tW^ provides a much more efficient trigger in terms of the W~ accompanying the top quark. 
So, the analysis presented in this section should be repeated with a more complete code, including 
all three channels. However, we think that in estimating the various efficiencies, the current version 
of PYTHIA is adequate. The distributions in the scaled energy variable Xk and in the transverse 
momentum pT of the K^s produced in the process pp — t- t/tX, and the subsequent decays t ^- W s 
and t ^ W b are shown in Fig. U] (upper two frames) for ^/s = 14 TeV. The corresponding 
distributions for the As from the decays t — t- VFs(— t- AX) and t — t- Wb{-^ AX) are shown in the 
lower two frames in Fig [H The scaled charged lepton energy distribution X£ = Eg/E^ct and the 
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FIG. 5: pp — )• t/tX at y/s = 14 TeV. Upper frame: scaled-i?^-energy distributions, dN/dxg, from t -^ W s(— >■ 
£^X) (dashed histogram) and t -^ W b{-^ i^X) (sohd histogram). Middle frame: Transverse momentum 
distributions of the is measured w.r.t. beam axis, dN/dpTi (i^ GeV), in the same production and decay 
processes as in the upper frame. Lower frame: Secondary decay vertex distributions in the variable r (in 
millimetres) for the two decay chains t -^W s (solid histogram) and t ^^W h (dashed histogram) , obtained 
by smearing the decay length with a Gaussian having an r.m.s. value 2 millimetres. 
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FIG. 6: pp ^ t/tX at yi == 14 TeV. Left: The normalised BDTD response, calculated by using the TMVA 
(see text). The signal (dark shaded) from the decay t -^ W s and the background (light shaded with dotted 
lines) from the decay t ~> W b are clearly separated in this variable. Right: Background rejection vs. signal 
efficiency calculated from the BDT(D) response. The two MVA methods yield very similar results. 

lepton transverse momentum are shown in Fig. [5] (upper frame), showing the distributions from 
the t — )• Wb{-^ £^X) (dashed histograms) and t — )■ Ws{^ £^X) (solid histograms). Transverse 
momentum distributions of the £^s, measured w.r.t. the axis, are also shown in this figure (middle 
frame). Finally, the secondary decay vertex distribution for the decays t — )• Wb (dashed histogram) 
and t — )• Ws (solid histogram) are shown in the lower frame in this figure. In plotting these 
distributions, we have smeared them, as before, with a Gaussian with a r.m.s. value of 2 millimetres 
and have taken into account the finite lifetime of the 6-quark, as stated in the case of the analysis for 
the process pp — t- tiX. These distributions provide an excellent discrimination between the signal 
t ^ W s and the dominant background t ^ W b events for the single top production process 
PP ~^ t/tX, qualitatively very much along the same lines as discussed earlier for the tt production 
pp — )■ tiX. As already stated, this information is used to build up a decision tree . In Fig. [6] (left 
frame), we show that the BDTD response function is very different for the signal (t — )• Ws) and the 
background (t — )• Wb) events also for the single top (antitop) production process. The background 
rejection vs. signal efficiency for the pp — )■ t/tX events is shown in Fig. [6] (right frame). The results 
for the signal efficiency and purity are shown numerically in Table II. The entries in this table show 
that also in single top production process a background rejection of 10^ can be achieved at a signal 
efficiency of about 5% to 7% to reach the SM-sensitivity of the CKM matrix element \Vts\- Thus 
we would expect 0.053 x 1.7 x 10^^ x 10^ = 90 signal events with a background of lO^'^ x 10^ = 10^ 
events giving a significance of 90/VToOO i.e. about 3a. 
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TABLE II: Tagging efficiencies (in %) for the single top production process pp — ?► t/tX, followed by the 
decay t — > Ws (signal) and t —> Wb{ background), calculated for an acceptance of 0.1% for the background 
at three LHC centre-of-mass energies. Two Gaussian vertex smearing (having an r.m.s. values 2 nun and 1 
mm) are assumed for calculating the displaced vertex distributions dN /dr. 



CTo 


7TeV 


lOTeV 


14TeV 


2 mm 


7.1 


6.5 


5.3 


1 mm 


21.7 


22.4 


19.9 



IV. SUMMARY AND OUTLOOK 

We have presented a case here to measure the matrix element \Vts\ from the top quark decays 
t — )• W^ s and its charge conjugate i — )• W~s, making use of the characteristic differences in the 
h- and s-jet profiles. We have concentrated on the VO {K^ and A) energy-momentum profiles 
emanating from the signal (t — )• Ws) and the dominant background {t — )• bW). This information 
is combined with the secondary vertex distributions, anticipated from the decays (5 — )• c — )• s), and 
the absence of energetic charged leptons in s-quark jets. An important parameter is the vertex 
resolution, for which we have used two values, (T(vertex) = 1 mm and 2 mm, assuming a Gaussian 
distribution. With these distributions, we train boosted decision tree classifiers, BDT and BDTD, 
and use the BDT(D)-response functions for the signal (t — t- Ws) and background (t — t- Wh) events. 
This information is used to study the background rejection versus the signal efficiency, which 
would enable to achieve typically 10% signal efficiency and a background rejection of 10^. Detailed 
studies are done at three representative values of the LHC centre-of-mass energies, ^/s = 7 TeV, 10 
TeV and 14 TeV. As the principal results (BDT(D) response functions and background rejection 
versus signal efficiencies) are very similar for all three energies, we present detailed results only for 
^ = 14 TeV. 

In this exploratory study, we have made some simplifying assumptions. In particular, we have 
used PYTHIA to undertake our analysis. The cross sections for the top pair [pp — )• tiX) and 
single top production {pp — )• t/iX) in PYTHIA can be adjusted to correspond to the theoretical 
precision currently available. However, the distributions and topologies, in particular for the single 
top (anti-top) production processes, will have to be correctly incorporated in a realistic simulation. 
Likewise, we have not attempted to define the s- and 6-quark jets using a modern jet algorithm. 
No attempt has been made at improving the training process by adding some more variables, like 



the b-jet shapes [36! ], which are known to have some discriminating power. We recall some of the 



important sources of uncertainties in our analysis: (i) predicted rates of the top quark production, 
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(ii) histogram shapes, (hi) integrated luminosities, (iv), efficiencies of the b- and s-quark tagging, 
reflecting in our study the relative efficiencies given in Tables I and II, and the uncertainty in 
rrit, though this effect mainly the CKM matrix element determinations form the single top (anti- 
top) production processes. All these experimental and theoretical refinements will have to be 
incorporated in the analysis of the LHC data to draw quantitative conclusions. In particular, 
background processes, most notably W + jets, Z + jets will have to be considered. Nevertheless, 
we have shown, in the ffist study of its kind, that a direct measurement of \Vts\ in top quark decays 
is, in principle, feasible at the LHC. The simulations presented here for 14 TeV correspond to an 
integrated luminosity of 10 fb~^. 

Alternative methods of determining the matrix elements \Vtd\, {Vtsl and \Vtb\ at the LHC are 
based on the single top (or anti-top) production at the LHC. One attempts to determine these 
inatrix elements from the cross section measurement by a simultaneous fit. This cross section is 
parametrised as a{pp — >• tX) = AdlVtdl"^ + AslVtsl"^ + Ab|Vjf,p (and likewise for a{pp — )• iX)), one 
then solves the cross-section for the CKM matrix elements, given the dynamical quantities A^, As 
and Afy. They, in turn, depend on estimates of the various electroweak processes in the single-top 
(or anti-top) production and on the parton distribution functions (PDFs). Typical estimates of 
the reduced cross sections at the LHC {^/s = UTeV) are: Ad = 766(253) pb, A^ = 277(172) pb, 
and At = 150(87) pb [37], where the numbers in the parenthesis refer to the production of anti-top 
at the LHC. Based on these estimates, one expects at the LHC As/Ab ~ 2 and Adj A^ ~ 5. These 
ratios depend on QCD and hence will not change if the weak interactions in the SM are modified 
by new physics. In the SM, one expects |Ftsp/|Vtbp ~ 1.6 x 10"^ and |FtdP/|^tbP ~ 6 x 10"^. 
In the example of realistic beyond-the SM physics that we are using to motivate these studies, 
these CKM matrix element ratios could be larger by a factor 4. We conclude that both in the 
SM, and in the four generation extension of it, the cross sections g(j)p — )• tX^ and a{pp — )• tX) 
are completely dominated by the ^6|Vf;,p term. Hence, this proposal does not have the desired 
sensitivity to measure the matrix elements |Vfd| and 1 14^1 at the level of theoretical interest. 



It has been recently suggested in [381] that one may improve the sensitivity to jVirf], if one refines 
the experimental analysis using the top quark rapidity distributions, which are different for the 
valence d-quark initiated processes as opposed to the sea 6-quark initiated processes [3^]. While 
of some value in exploring iVj^j, still the sensitivity of this method is far from the expected value 
of \yid\ by an order of magnitude. Moreover, as the s-quark and the 6-quark are both sea-quarks 
in the proton, the top quark rapidity distributions do not provide an improved determination of 
\yts\ from the single top production process. Hence, our method based on the top quark decay 
characteristics to determine \yts\ complements the existing proposal. 

Finally, we remark that the ratio of the CKM matrix elements (] 14(^1^ -|- ] Visl^)/lVf;,]^, that can 
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be obtained by measuring the ratio Rt, defined in the introduction, through the number of events 
with zero-, one-, and two 5-tags in the process pp — )■ tiX, can be combined with the determination 
of the ratio |VisP/|Vjf,p discussed here, to constrain (or measure) the quantity iV'tdP/l^tbP- 
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